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Abstract. Beliakova gave in [2] a method for computing link Floer homology. We extend her 
approach to Z coefficients and implement it in a publicly available computer program. We review 
the main programming and optimization methods used. Our program is then used to check that 
the Floer homology of a prime non-alternating knot with < 11 crossings has no torsion. 



1. Introduction 

Knot Floer homology HL was introduced in [12] by Ozsvath and Szabo and independently in 
|15| by Rasmussen. In this article, we will restrict ourselves to the study of HL, a simplified version 
of HL. The knot Floer homology HL(K) of a knot K is the homology of a bigraded complex (C, d) 
with Maslov grading m and Alexander grading a. 

HL(K) = HL m (K, a) 

We will work with either Z or Z/2Z as coefficient ring. Knot Floer homology, even in its simplest 
version HL{K) over Z/2Z, gives much information about a knot: The Seifert genus g{K) of K is 

arg max( HL{K, a) ^ 0) 

a 

(see [H]). The knot K is fibered if and only if iankHL(K,g(K)) = 1 (see [5] and [9]). It is also 
possible to extract from HL a bound for the slice genus of the knot, see |13j . 

The first combinatorial way for calculating knot Floer homology was given in [10] for Z/2Z 
coefficients and was then extended for Z coefficients in [8]. Those papers explain how to construct 
a complex (C, d) from the grid diagram (also called rectangular diagram or arc presentation) of 
K, so that HL(K) = H*(C) <g> F® (n_1) (see Theorem [23] for the definition of the free module 
V). However, the number of generators in C increases very fast with the size of the grid diagram, 
making the practical computation of this complex difficult. A computer program extracting knot 
Floer homology from this complex was written by Baldwin and Gillam (see [1]). This program 
manages to alleviate the problem of the number of generators by taking advantage of the sparseness 
of the matrix describing the boundary map d. 

In [2], Beliakova proposed a new complex (C\ ong , d\ ong ) for knot Floer homology, which we will 
call the long oval complex. She explained how it is homotopic to a complex (C s hort> cohort); which 
we will call short oval complex. The main interest of the construction is that the short oval complex 
usually has a much smaller number of generators than the MOS complex. Those two complexes 
are described in Section 2. 

In Section 3, we extend Beliakova's complexes, originally defined with Z/2Z coefficients, to Z 
coefficients. In Section 4, the formula used to calculate the Alexander grading of generators in oval 
complexes is proved. 

Section 5 explains how our program computes the short oval complex and extracts its homology. 
It also describes algorithms and optimizations that could have a larger interest. Many could, for 
example, be used to work with the MOS complex. 
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1.1. Main results. We developed a program that can be used to determine the HL with Z or 
7LjTL coefficients of knots with less than 13 crossings efficiently and can be used to determine 
fiberedness and the Seifert genus of even bigger knots... Using our program, we show that the HL 
of prime non-alternating knots with less than 12 crossings contains no torsion and checked the HL 
computations given in pQ. 

1.2. Acknowledgements. I would like, above all, to express my gratitude to Anna Beliakova for 
the original impetus of the project, numerous interesting discussions, advice and ideas. I also would 
like to thank Dror Bar-Natan for good advice and his support in making my program a part of the 
Knot Atlas. The resources of the Knot Atla^ and of the knot data tables of Alexander Stoimeno^ 
were very useful to me. Some of the figures in this article are courtesy of Anna Beliakova. 

2. Definitions of the complexes 

2.1. MOS complex. A grid diagram (or grid diagram on the plane) of complexity n is a n x n 
square grid on the plane, with each square decorated by an X, an O or nothing, such that each 
column or row of the grid contains exactly one X and one O. The set of all O (respectively all X) 
is called O (respectively X). (See Figure [TJ) 



r-Xr-r--- 



6* 



-r-V-r 



-Kj 



Figure 1. A grid diagram for the 3i knot. The white points represent 
black points O. 



and the 



We call the center of a decorated square a puncture. More information about grid diagrams can 
be found in [3J. 

A grid diagram represents a link in the following way. If we draw a line segment between pairs 
constituted of an X and an O sitting in the same row or in the same column, we obtain a figure 
in the plane. Assuming that vertical segments always pass over horizontal segments when they 
intersect, we can see this figure as a projection of a link. We say that the grid diagram represents 
the same link as this projection. It is not difficult to see that a knot can always be represented by 
many different grid diagrams. 

Starting with a grid diagram on the plan, by identifying the uppermost and bottommost lines 
of the grid and the leftmost and rightmost lines, we obtain a decoration of a toroidal grid called a 
grid diagram on the torus. After this operation, the vertical and horizontal lines of the grid become 
meridional and longitudinal circles. Each grid diagram on the torus of complexity n is associated 
in a natural way to n 2 grid diagrams in the plane. 



^Our program is accessible on our homepage "http:/ /www. math. unizh.ch/assistenten/jdroz" or via the Knot Atlas 
'http: /katlas. math. toronto.edu/wiki/Main_Page" 
2 http: /katlas. math.toronto.edu/wiki/Main_Page 
3 http:/ /www. kurims.kyoto-u.ac.jp/ stoimeno/ptab/index.html 
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Definition 2.1. Let D be a grid diagram on the torus of complexity n, and D' one of its analog 
in the plane. The MOS complex (C(D),d(D)) is defined in the following way: C{D) is the vector 
space over Z/2Z, which has as basis the unordered ra-tuples of intersection points of circles of the 
grid that contains exactly one intersection point on every meridional or longitudinal circle. The 
boundary map is defined as follows: 

j/6Rocto(x) 

We define Recto (x) to be the set of generators y such that: 

• The generators x and y have exactly n — 2 intersections in common. 

• The four intersection points that are either only in x or only in y form the four corners of 
a rectangle R. In D' , the upper right corner of R is in x. 

• The rectangle R is empty in the sense that it doesn't contain a puncture of the grid diagram 
or any intersection point of x. 

We often think of the boundary map as a matrix with entries indexed by the generators of the 
complex: 

d{x) = ^2 d *,y ' V- 
y 

In other words, d = (d xy ) for x, y generators. 

Definition 2.2. let S,T C M 2 be two finite sets of points in the plane. I(S, T) is defined as 
the number of pairs (a, b) G S, (c, d) G T with a < c and b < d. We also define J(S, T) = 
(I(S,T) + I(T, S))/2. We consider the natural bilinear extensions of / and J to formal sums of 
intersection points. 

Definition 2.3. A generator g € C has Maslov grading M(x) = I(x, x) — I(x, O) — 1(0, x) + 
1(0, 0) + 1. It will also sometimes be practical to speak of the Maslov grading of an arbitrary set 
of points. 

Definition 2.4. A generator g 6 C{D) for a grid diagram D of complexity n has Alexander grading 
A(x) = J(x - (O + X)/2, X - O) - (n - 1) /2. 

It is also possible to write the Alexander grading of a generator x in a diagram D as the sum of 
the winding number of D around the intersection points pdi and a constant depending only on 
D. Let a(p) denote the average of the winding number of the knot projection represented by D in 
a small ball around p. 

Despite a small difference in the definition of a(p), the formula above is equivalent to the formula 
for the Alexander grading given in [10j . 

Theorem 2.5. (C.Manolescu, P.Ozsvdth, S.Sarkar) Let D be a grid diagram of complexity n of 
the knot K, (C(D),d(D)) the MOS complex constructed from D. 

• The chain complex (C(D),d(D)) is bigraded by the Alexander and Maslov gradings. The 
boundary map preserves the Alexander grading and decreases the Maslov grading by 1. 

• The homology of H*(C(D),d) is HL(K) (g> y®*-™ -1 ^ where V is a vector space with basis 
composed of one vector of Alexander and Maslov gradings -1 and one vector with Maslov 
and Alexander gradings zero. 
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Figure 2. A grid diagram for the 52 knot with a collection of long ovals. 

2.2. Beliakova's complexes. We will now construct two complexes that have the same homology 
as the MOS complex, and for which a theorem almost identical to Theorem 12.51 is proven in [2]. 
The construction of the long oval complex (C\ ong (D), d\ on g(D)) also starts with the grid diagram 
D{K) of the knot K. But this time, the construction is done in the plane. Let us call n the 
complexity of the grid diagram D. We begin by drawing n — 1 long thin vertical (respectively 
horizontal) ovals around all but one pair of punctures with identical first (respectively second) 
coordinates (see Figure [2]). We call the vertical ovals a = {a\, . . . , ct n -l} an d the horizontal ovals 
(3 = . . . , fl n -i}. While drawing those ovals, we also require that at least one puncture is in the 
unbounded component of M 2 /(ai U ...a n _i U fi\ U .../3 n _i). 

Definition 2.6. Let D be a grid diagram of complexity n in the plane. The vector space over 
Z/2Z called C\ ong {D) has as basis elements the sets of n — 1 intersection points between vertical and 
horizontal ovals with exactly one intersection point on each oval. An Alexander grading on those 
generators is defined by exactly the same formula as in the case of the MOS complex. The Maslov 
grading M(x) of a generator x is defined to be I(x,x) — I(x,Q) — I(Q,x) + 1(0,0). (The same 
formula as for the MOS complex minus one). The boundary map d\ ong (D) is defined as follows: 

yeRecto(a;)U?/eBigonQ(a;) 

The set Recto (x) contains the generators y such that: 

• The generators x and y have exactly n — 3 intersection points in common. 

• The four intersection points that are either only in x or only in y form the four corners of 
a rectangle R. The sides of R are arcs of ovals. The upper right corner of R is in x. 

• The rectangle R is empty in the sense that it doesn't contain any puncture of the grid 
diagram or any intersection point of x. 

The set Bigon (a:) contains the generators y such that: 

• The generators x and y have exactly n — 2 intersection points in common. 

• The two intersection points that are either only in x or only in y constitutes the two corners 
of a bigon L. The sides of L are one arc of horizontal oval and one arc of vertical oval. 
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Figure 3. A grid diagram for the 52 knot with a collection of short ovals. The 
blue points represent a generator. 

A counterclockwise rotation around L, along the arc of the horizontal oval, leads from the 
corner in x to the corner in y. 
• The bigon L is empty in the sense that it doesn't contain any puncture of the grid diagram. 

The long oval complex can be reduced, by a sequence of homotopies (hi, hi), to a much 
smaller (in terms of the rank of Ci on g(-D)) complex, which we will call the short oval complex 
(Cshort(-D)j cohort (see Figured]). This sequence of homotopies corresponds to a progressive 
shortening of the curves a and the curves /3. This shortening is an homotopy written (a,(3)t for 
t € [0, 1 + 1] such that there is no t for which a curve meets a puncture. Among such homotopies, 
we choose one that minimizes the number of intersection points of the curves at t = I + 1. We 
choose this progressive shortening such that one pair of intersection points disappears for each 
t G {1,...,/} (see Figure [J] for an example of the corners of a bigon disappearing). To each (a,(3)t 
t E [0,1 + 1], we associate a complex (Cm (D), <9^j (D)). The complex Ck(D) is generated by sets 
of n — 1 intersection points of and (3 k such that each intersection point is on exactly one oval. 
The boundary maps dk(D) are defined inductively. 




Figure 4. The disappearance of a bigon during an homotopy from the long oval 
complex to the short oval complex. 
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Definition 2.7. Let p± and p2 be the intersection points that disappear at time t = k, assume 
{pi} has a bigger Maslov grading than {p2}- For x a generator of Ct{D), let rj(x) be either 0, if x 
does not contain P2, or, if x does contain P2, a generator identical to x except that it contains p\ 
instead of P2 ■ 

d k = 7T O (d k -l + <9fc-l O 7] O o 1 

Where tt is the natural projection from Cfc_i(-D) to Cfc(D) that sends generators containing pi or 
P2 to zero and 1 is the natural injection from Ct{D) to Ck-i(D). 

The homotopy /i/t is an homotopy between (Ck-i(D), dk-i(D)) and (Ck(D), dk(D)). It corre- 
sponds to the "cancellation" of two intersection points. The existence of these homotopies is Lemma 
2.1 of [2]. 

It is noteworthy that this construction of (C s hort(-D)> cohort is not canonical. The differential 
cohort depends on the whole (a,f3) t (for t £ [0,1 + 1]) and not just on (a,(3)i + i. 

Theorem 2.8. (Beliakova) LetD be a grid diagram of complexity n of the knot K , (C s hort(-D); dshorx,{D)) 
a short oval complex constructed from D. 

• The chain complex (C s hort(-D)i <9short(-D)) is bigraded by the Alexander and Maslov gradings. 
The boundary map preserves the Alexander grading and decreases the Maslov grading by 1. 

• The homology H*(C short (D),d shovt (D)) is HL(K) ®V^ n ~ l \ where V is a vector space 
with basis composed of one vector of Alexander and Maslov gradings -1 and one vector with 
Maslov and Alexander gradings zero. 

3. Sign Assignment 

We start by defining an extension (C{ ong (D), d[ (D)) of the long oval complex, the homology 

of which is HL over Z. The complex (C[ (D), d' long (D)) will be said to be a sign assignment over 
the long oval complex. 

Definition 3.1. Let (C',d') be a complex with Z as coefficient ring and basis b' = {b' l7 . . . ,b' n }. 
Let (C,d) be a complex with Z/2Z as coefficient ring and basis b = {61, . . . ,b n }. Let us assume 
that the matrices {d[ j) and (dij) represent the boundary maps in the basis b' and b. We say 
that (C, d') is a sign assignment on (C, d) if the following two conditions are true for all integers 
1 < i, j < n. 

. d l)3 = & id = 

• dij ^ => d'i j = +1 

Our approach is completely analogous to the one taken in [8] for the same purpose. We assume 
for the following definitions that the curves in afl/3 are ordered. We say that a point sits positively 
on an oval if it sits on the right side of an a curve or the upper side of a f3 curve. 

Definition 3.2. Let x and y be generators of the long oval complex. We define Sign(x,y) in the 
following way: 

• If y £ Recto(x) (we say that x and y are connected by a rectangle): Let (0,6) (respectively 
(c, d)) denote the coordinates of the leftmost (respectively rightmost) point that is in x but 
not in y. Let D(x,y) be the number of points p = (pi,P2) in % such that a < p\ < c and 
P2 < b. In other words, D(x, y) is the number of points below the rectangle between x and 

y- 

Sign(x y) = ( — iy( x '{( x l> x 2)£x\x2<d})+D(x,y)-(I(x,{(x 1 ,x 2 )£x\b<x 2 <d})+l) 

• If y G Bigon (x) (we say that x and y are connected by a bigon): Let E be the oval on 
which x and y have intersection points on opposite sides (the intersection points are placed 
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symmetrically with respect to the great axis of E). Let pre(x,y) be the number of ovals 
that come before E in the ordering and on which an intersection point of x sits positively. 

Sign(x,y) = (_l) / (*.*)+P«(*.») 

• Otherwise, Sign(x, y) is undefined. 

Definition 3.3. The complex C[ Q (D) is a free module with Z coefficients on the same generators 
as Ci ong (£>). The boundary map df ong is the linear map equal to Ej /G Rect (x)uBigon (x) Sign(x,y) • y 
for each generator x of C. 

For generators (x,y) between which there is a rectangle (y G Recto(x)), our formula is taken 
from 0. 

Lemma 3.4. The boundary map d[ verifies the condition d'f ong = and thus (C[ ong (D), d[ (D)) 
is a chain complex. 

Proof. Let x and z be two generators of C[ Q (D) with Mix) — 2 = M(z). The number of generators 
y such that y G Recto(:c) U Bigon (a;) and 2 6 Recto (y) U Bigon (y) is either or 2. If those two 
generators exist, we call them y\ and 2/2 • The lemma is equivalent to the claim that when y\ and 
2/2 exist, Sign(x, yl) • Sign(x, y2) ■ Sign(yl, z) ■ Sign(y2, z) = —1. We must check this in three cases: 

• First case: 2/1,2/2 G Recto (sc). This implies that z G Recto (2/1) and 2 G Recto (2/2)- Therefore, 
the claim follows from [8] Section 4.1, where an analogous claim for the MOS complex is 
proven. 

• Second case: y\ G Recto(:r) and 2/2 G Bigon (x). This implies z G Recto(2/2) and z G 
Bigon (yi). Since the formula for Sign for generators connected by a rectangle only depends 
on which pairs of ovals have common intersection points, we have Sign(a;, 2/1) = Sign(2/2, z). 
Since pre(x,2/2) = pre(j/i,^) and I(x,x) ^ I(yi,yi) (mod 2), Sign(x,2/2) = -Sign(2/i,z) 
and the claim follows. 

• Third case: 2/1,2/2 G Bigon (x). This implies z G Bigon (2/i) and z G Bigon (2/2). We 
assume w.l.o.g. that the oval containing the bigons connecting x and 2/1 and 2/2 and z 
comes before the oval containing the bigon connecting the other pairs. We then have 
pie(yi,z) pre(x,2/2) (mod 2) and pve(y2,z) = pie(x,y±). Because for a generator g, 
I(g, g) does not depend on which side of the ovals the intersection points of g are, we have 
I(x,x) = 1(2/2,2/2) = I(yi,yi)- The claim follows. 

□ 

Theorem 3.5. Let K be a knot with rectangular diagram D. The homology of (C[ ong (D), d[ Q (D)) 
is HL(K) <g> T/®( n_1 ) (see Theorem\2J& for the definition ofV). 

Proof, (sketch) We begin by extending our sign assignment for (C\ ori g(D), c\ ong (D)) to a sign as- 
signment for a modification of the long oval complex called (Cf u n, <9f u n) that has the full HL(K) <g> 
V<a(n-i) ( and nQt hl(K)) ® y®^" 1 )) as homology. In [8] Section 2.3, the complex for HL(K) ® 
y®{n-i) ig » extractec i" from the more complicated complex for HL in the same way that the long 
oval complex can be extracted from (Cf u u, 9f u n). 

We then prove that there is, up to quasi-isomorphism, only one sign assignment (Cf ull , d' iun ) 
with dLy[ = 0. This is analogous with [8] Section 4.1. The existence of a sign assignment for 
(Cfull, dfuii) with homology HL(K) ® V^ n ~^ is a consequence of the analytical theory of link Floer 
homology (see [12] for the general theory and [2] for the special case of the long oval complex). 
Since our (extended) sign assignment verifies the condition dL u = (this is proved in the same 
way as Lemma I3.4H . it must be quasi-isomorphic to the one coming from the analytical theory. 
Quasi-isomorphic complexes have identical homologies. The theorem follows by restricting the 
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sign assignment on df n \\ to a sign assignment on <9iong and noticing that because the homology of 
(Cfuli; d U) is HLiK^V^-V, the homology of (C{ ong (D), d[ ong (D)) will be HLiK^V®^. □ 

The short oval complex over Z/2Z can now be extended to a complex with Z coefficients. We 
get a "signed" short oval complex by applying to the new long oval complex a tiny modification of 
the inductive construction we used in the modulo 2 case. Instead of Definition 12.71 we take: 

Definition 3.6. Let p± and p2 be the intersection points that disappear at time t = k, assume 
{pi} has a bigger Maslov grading than {p2}- For x a generator of Ck(D), let n(x) be either 0, if x 
does not contain P2, or, if x does contain P2, a generator identical to x except that it contains p\ 
instead of P2 ■ 

d k = 7ro (<9 fe _i - <9 fc _i ot]o <9 fc _i) o 1 
Where tt is the natural projection from Ck-i(D) to Ck(D) that sends generators containing p± or 
P2 to zero and 1 is the natural injection from Ck{D) to Ck-i(D). 

And instead of lemma 2.1 of [2] we use: 

Lemma 3.7. The complexes Ck{D) and Ck-i(D) of Definition \3.b] are homotopy equivalent. 

Proof. This lemma can be proved either by "adding" signs to the proof in [2] or by noticing that the 
generators containing a corner of the bigon can be matched in the sense of discrete Morse theory 
and that Ck(D) is Ck-i(D) with the matching collapsed (see [6]). □ 

Our new short oval complex is not a sign assignment over the one defined in the previous section 
since some coefficients of the boundary map can be of absolute value greater than one. 

4. The Alexander grading 

In this section, we give a proof of Formula [231 for the Alexander grading of the long oval complex. 
Let us first formulate a technical lemma. 

Lemma 4.1. The minimal number I such that a permutation a is a product of I transpositions is 
written sign(cr). Let M be a matrix of size n with all entries in its first row and first column equal 
to 1. Then, (T7]j holds. 

(1) det(M) = det(Mij) 

l<i ,j<n 

Where Mij denotes the minor of M obtained by removing the row number j and the column 
number i. 

Proof. Let us denote the entries of M by rrnj for 1 < i,j < n. Once we apply 

(2) det(A) = J] a lMl) ...a nAn) (-lf^ 

cr€S n 

to both sides of dU), we get an equality between polynomials. 

(3) E m M(D- m n,a(n)(-l) Sign(<T) = 

(4) J] (-1)^ X ( _l)si g n(<x) + tt {fc :0<( fc -i)( CT(fe )-i)} TJ mka(k) 
Ki,j<n a£'S n ,a(i)=j l<k<n,k^ii 

We will prove it by comparing coefficients of monomials on both sides, while setting mij = 1 and 
m^i = 1. We classify the monomials appearing in ([3]) in two types: ?7ii j(T (i)--- m n,o-(n)( — l) slgn( - <T ) for 
a G S n with er(l) = 1 and "^1^(1) ■■■ m n,a{n)i~ l) slgn ^ with cr(l) 7^ 1. The sum of the monomials of 
the first type are exactly the monomials appearing in (|4|) coming from det(Mi j i). The monomials 



EFFECTIVE COMPUTATION OF KNOT FLOER HOMOLOGY 



9 



of the second type all appear in exactly three times, two times with the same sign as in ((3]), 
one time with opposite sign: (we use G to denote that a monomial is a summand of a polynomial 
written in canonical form) 

(5) m 2iCT{2) ...m niCT(n) (-l) si ^ G det(M ljff(1) ), det(M a -i (1)jl ), - det(M CT - 1(1)/j(1) ). 
All monomials of ((3|) are now accounted for in (j4j). But monomials of the type: 

(6) ± mi,a(i) ■■■m i _ l ^ i _ 1) m i+l ^ i+1) ...m n ^ n) G det(M i a ^) 

for some i ^ 1 with a(l) ^ 1 and a(i) ^ 1, remain in Setting W to be equal to a composed 
with the transposition exchanging i and <t -1 (1), we notice that 

(7) (_ 1 )sign(a)+tt{fcO<(fe-i)( e T(fc)-a(i))}+i+a(i) "Q mfe (T(fe) ( G ± det(M i)Cr(i) )) 

l<k<n,k^i 

is equal to 

( 8 ) ^_ 1 ^ ) l+sign(a)+t){fc:0<(fc-a- 1 (l))(a(fe)- ( T(i))}+ CT - 1 (l)+<T(i) m kw{k) 

( G ±det(M CT -i (lWi) )). 

Therefore, those monomials cancel each other too. □ 

Theorem 4.2. The Alexander grading on generators of the long oval complex on a grid diagram 
D is given by the same formula as the Alexander grading on generators of the MOS complex: 



Here, n is the complexity of D, a(p) the winding number of D around p and O one of the two types 
of punctures. 

Proof. It follows from the analytical theory of knot Floer homology, that the relative Alexander 
grading on pairs of generators must be defined in the same way in the two complexes (see |10| and 
[2] ) . We prove the theorem by showing that the Euler characteristics of the long oval complex and 
the MOS complex are equal: 

X(Clong) = X(CMOS)- 

Since the Alexander polynomial is never null, this fixes the absolute Alexander grading. 

Let M be an n x n matrix, with entries m^j equal to t power the winding number of the knot 
around the grid diagram intersection point of coordinates (i,j). The Euler characteristic of the 
MOS complex can be calculated by taking the determinant M and multiplying by a constant f{D) 
depending only on D. 

Drawing the diagram that gives rise to the long oval complex, as in Figure [5j we can associate 
sets of intersection points of the grid with oval complex generators. Working on the torus, we do 
this, as in Figure EJ by replacing intersection point of ovals by the nearest intersection point of the 
grid. It is easily checked (by pairing generators of opposite contributions, see Figure [6]) that the 
total contribution of the generators associated with sets having two points with the same vertical or 
horizontal coordinate is zero. But each set of intersection points of the grid with one point at most 
on each vertical or horizontal line of the grid is the associate of exactly one long oval generator. 

Let S be any of those sets of n — 1 points and g its associated generator. The set S misses exactly 
one row j and one column i of the grid diagram. The set S can be seen as a monomial of the minor 
Mi j, the point of coordinate (k, I) in S corresponding to the entry m kt i of M. The Alexander index 
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Figure 5. A generator on a set of long ovals and the set of nearest intersection 
points on the grid. 
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Figure 6. A set of n — 1 intersection points with two points on the same column is 
associated here with two generators of the long oval complex. Since those generators 
have the same Alexander grading and since their Maslov gradings differ by one, they 
cancel each other in the Euler characteristic. 



of g is exactly the value of the monomial of the minor Mj j. Moreover, the sign of the monomial in 
(— ■ Mij is the Maslov index of g. Collecting those facts, we get: 

x(ao„ g ) = (E(- i r-Ei^n a g i^)-/(^)= E (-ly^'det^)-/^ 
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Using Lemma [H 

= det(M) • f(D) = x(Cmos). 

□ 

5. The implementation 

For reasons of simplicity, we describe our program only for Z/2Z coefficients. Passing to Z 
coefficient doesn't create any new problem. Our program computes the knot Floer homology 
HL(K) of a knot K by constructing a short oval complex (C s i lort (D(K)), d s i 10Ti (D(K)) and taking 
its homology. Our program therefore contains three main parts: 

• A function taking as input a knot (in braid representation usually) and giving as output a 
rectangular diagram with a set of ovals. The rectangular diagram and the ovals are chosen 
in order to minimize the number of generators of the short oval complex constructed from 
the ovals. 

• A function listing the generators of the complex and calculating their Alexander and Maslov 
gradings. 

• A function that computes the boundary map. Since most of the running time of the program 
is spent by this function, it is heavily optimized. 

• A function calculating the homology of the complex. This function gives in fact HL(K) (g) 
y®(n-i)_ Extracting HL(K) from it is easy. 

The first three parts are explained in the next three subsections. The greatest task in the computa- 
tion is to calculate the differential cohort- This differential, being a linear map, could be represented 
as a huge matrix Mq. However, because the differential preserves the Alexander grading and de- 
creases the Maslov grading by one, the differential can be represented by a set of matrices of much 
more reasonable size, one for every pair of Alexander and Maslov gradings. 

Our short oval complex is a direct sum of complexes with uniform Alexander gradings. 

Cshort = C s hort 
(i 

Here, C° hort is the subcomplex of C s h or t generated by generators of Alexander grading a. Since we 
compute an homology of the form HL(K) <g> y®( n_1 ) with V a vector space spanned by two vectors 
of different Alexander gradings, we don't even have to compute the whole complex Cohort- F° r an y 
set A of n— 1 Alexander gradings, the complexes C" hort can be deduced from the rest ((J)^^ C" hort ). 
For the Alexander gradings we choose to ignore, we don't have to compute the generators or the 
boundary map. Naturally, we choose to ignore the Alexander gradings that would be the most 
difficult to compute: those with many generators. Since the numbers of generators of the possible 
Alexander gradings is very unevenly distributed, a lot of time is saved. 

5.f . Generating a grid diagram. Generating a grid diagram representing a knot, either by hand 
or automatically, is not difficult (see Figure [7] for an example). However, the grid diagram obtained 
usually has much higher complexity than what is necessary. Therefore, our program simplifies 
knot diagrams before using them to construct complexes. This is done by exploring the set of 
equivalent rectangular diagrams during a fixed amount of time and picking the diagram of least 
complexity generated during this period. A set of simple moves (called Cromwell moves, see [3] or 
[1] ) analogous to Reidemeister moves are used to this effect (see Figure [8]) . A cycling move is a 
circular permutation of the rows or column of a grid diagram. A stabilization move is the merging 
of two rows (or column) containing adjacent decorated squares accompanied by the deletion of the 
column (respectively the row) containing the two squares. A destabilization is simply the inverse 
of a stabilization. A castling move is an exchange of adjacent rows or columns of the grid diagram. 
(Castling moves are only possible if the decorations on the adjacent rows or columns are in a certain 
order.) The search is optimized in four ways: 
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FIGURE 7. A representation by braid closure (top left) of the figure eight knot is 
transformed in four simple steps in a grid diagram (bottom right). 



• By caching (memorizing) already generated diagrams and exploring only diagrams of min- 
imal complexity. This enables us to avoid exploring the same diagrams many times. 

• By considering as identical diagrams that can be obtained from one another by cycling 
moves. By reducing the space we have to explore, this makes caching much more beneficial. 
This is more or less equivalent to having the diagrams sitting on the torus instead of on the 
plane. 

• By exploring only diagrams of minimal complexity. In other words, moves increasing the 
complexity are not considered and each time a diagram of smaller complexity is discovered, 
the search begins anew from this diagram. This does not only gain time but also makes the 
memory consumption of the caching more reasonable. Moreover, searching in a monotonic 
way probably doesn't miss too many opportunities (see [4]). 

• By generalizing the stabilization move, so that simplifications can be discovered faster. 
For example, sequences of castlings in the same direction followed by a stabilization are 
considered as one operation. 

Once we have a minimized grid diagram D of complexity n, the next step is to embed this diagram 
in the plane and to choose 2n — 2 ovals on it. The generators of the complex (C s hort(-D)) cohort (-D)) 
that we want to calculate are sets of intersection point of ovals. In order to minimize the number 
of generators, we simply generate every possible embedding of the diagram and every choice of 
(shortened) ovals on them and pick the set of ovals that has the smallest number of intersection 
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Figure 8. The three Cromwell moves: (a) a cycling move, (b) a stabiliza- 
tion/destabilization move and (c) a castling move. 

points. Of course, having a smaller number of intersection points of ovals does not guaranty having 
a smaller number of generators, but counting intersection points is easier than counting generators. 

5.2. Listing generators. Although constructing the generators is not very difficult, some care 
must be taken to compute their Alexander and Maslov gradings efficiently. To compute the Alexan- 
der gradings quickly, we tabulate J({p} — (0 + X)/2,X — O) beforehand for each intersection point 
p. Then for a given generator g, the Alexander grading is simply a sum of n — 1 tabulated values 
and a precomputed constant. 

The case of the Maslov grading is more difficult. The naive algorithm for computing I(x, x) 
for a set x C M. 2 , takes time 0(|x| 2 ). By using a divide and conquer approach (dividing M 2 by 
horizontal and vertical lines that split x in two for example), we can get an algorithm running 
in time 0(\x\ log |x|). However, in this case, the constant is quite bad. Our solution consists in 
generating the generators in two stages. Each proto-generatorg^ of the first stage giving raise to 
many generators in the second stage. This is done in a way that enables the worst part of the 
calculation, I(x,x), to be done on the proto-generators and to be inherited by the generators. 
The rest of the calculation (— I(x, O) — 1(0, x) + /(O, O)) uses the same tabulation method as the 
Alexander grading. 



The proto-generators simply contains a pairing of the vertical ovals with the horizontal ones. They are therefore 
between 2™" 1 and 4™" 1 less numerous than the generators. 
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Since, the fiberedness and the genus of a knot only depend on the part of its homology with high 
Alexander grading, it is noteworthy that the set of generators S of Alexander grading > c can be 
listed in time 0(|5*| • n 3 ) (at least when \S\ » n). This comes from the fact that generators can 
be seen as matchings in a weighted bipartite graph between vertex representing the horizontal and 
vertical ovals. The weights can be chosen so that the Alexander grading of a generator is simply (up 
to a constant) the weight of the matching. Listing matchings with weights bigger than a constant 
is done by using a backtracking search that tests for cuts by the Hungarian algorithm for bipartite 
weighted matching (see [7j). 

5.3. The boundary map. Computing the boundary map means computing the entries d XjV of 
the matrix representing it as a function of the generators x and y. When d xy ^ 0, we sometimes 
say that x and y are connected by the boundary map. We do not have a single method that always 
computes d x ^ y very quickly. However, d x ^ y is zero most of the time and, when d X)V is zero, we can 
usually determine that quickly. We therefore use two methods to compute d x>y : 

• A fast method to check if there is a "domain" between two generators. The existence of a 
domain is a necessary condition for d X)V to be non-zero (see l5.3.11 below where the definition 
of domains is given). 

• A method that calculates d XjV , which is inefficient, but which our program only calls when 
there is a domain between x and y (see I5.3.2I) . 

5.3.1. Domains between generators. Let D be a grid diagram on the plane with sets of vertical 
and horizontal ovals a and (3. The ovals divide the plane in many bounded connected com- 
ponents and an unbounded component, which we call pieces. The bounded components are 
bigons, rectangles and triangles. Domains are multisets (formal positive integral linear combi- 
nations) of pieces. For a domain D and an intersection point p of the ovals, the multiplicities 
in D of the four pieces adjacent to p listed in counter clockwise order starting from the upper 
right will be written a\ (p) , 02 (p) , 03 (p) , 04 (p) . The corner index c(p) of an intersection point p is 
a i(p) + a 3(p) ~ a 2{p) — 03 (p). The sum of the pieces in an oval is called a periodic domain. 

Lemma 5.1. Given a set of transverse vertical and horizontal thin ovals in the plane, the only 
domains D such that every intersection point of the curves has corner index zero are the sums of 
periodic domains and a multiple of the unbounded piece. 

Proof. The result is obvious when the ovals have no intersection points. The lemma is then proved 
by induction on the number of intersection points. □ 

Theorem 5.2. Given a collection of ovals on a grid diagram, the associated complex (C,d) and 
two generators x,y £ C. If d X)V 7^ 0, there is a unique domain D with the following properties: 

• The corner index at intersection points contained in x but not in y is 1. 

• The corner index at intersection points contained in y but not in x is -I. 

• All other intersection points have corner index zero. 

• Pieces of D containing punctures have multiplicity zero in D. 

Proof. For long oval complexes, the existence and uniqueness of D follow easily from the definition 
of the boundary map. For short oval complex, existence is proved by induction on the sequence of 
homotopies. 

Uniqueness of domain between generators in short oval complexes follows from the four conditions 
on D. Since the unbounded piece contains a puncture, the fourth condition implies that the 
multiplicity in D of the unbounded piece is null. Therefore, using Lemma 15.11 and the four three 
conditions on D, we get uniqueness of D up to addition of linear combination of periodic domains. 
However, at least one puncture is in only one periodic domain po an d all the periodic domains can 
be arranged in a sequence po, . . . ,P2n-2 such that pi and pi + \ have one puncture in common. Using 
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the fourth condition on D, this sequence enables us to show by induction on i that the multiplicity 
of pi in the difference of two domains verifying the four conditions must be zero. □ 

To determine the existence of a domain is therefore equivalent to solving in integers a system of 
linear equations and inequalities, the unknown being the multiplicities of the pieces in D. Luckily, 
the equations coming from the conditions on D in Theorem 15.21 always have a unique solution, 
even over Q. We can therefore decide efficiently if there is a domain between generators, by first 
solving a system of linear equations over Q and secondly checking that the solution obtained is 
non-negative and integral. 

5.3.2. Computing d Xj y. The boundary map d (for both long or short oval complexes) can be seen 
as a graph with vertices representing generators and edges representing the non-zero entries of the 
matrix d. Given an homotopy between long ovals and short ovals, the graph representing the short 
oval complex can be constructed by a sequence of modifications on the graph representing the long 
oval complex. The arguments used in this section are typical of algebraic Morse theory [6]. 

Each time a pair of intersection points forming the corners of a bigon disappears (see Figure H]), 
the graph is modified in the following way: 

• All generators containing one of the disappearing intersection points are deleted. 

• Let a and b be two generators connected by the bigon that disappears during the homotopy: 

b € Bigon (a). 

For every generators x and y with d X: b = 1 and d a>y = 1, we put one edge between x and y 

exactly if it wasn't there before. 
This is a simple translation of Definition 12.71 in graph theoretic language. Intuitively, it means that 
the edges of the graph of a short oval complex correspond to paths in the long oval complex graph. 

Thinking in terms of graphs enables us to give a clearer condition for d xy = 1 (mod 2) in the 
short oval complex (see Figure [9] for an illustration of the effect of an homotopy on the graph 
representing the boundary map). 

Theorem 5.3. Let (C\ ong , c\ ong ) be a long oval complex. Let (C s hort> cohort) be a short oval complex 
obtained by an homotopy (a,f3)t for t G [0, Z + 1] . Let G be a graph representing {C\ ong , c\ ong ). We 
associate to each edge of G, between generators connected by a bigon that disappears during the 
homotopy, the time at which it disappears. Then x and y are connected by the boundary map of the 
short oval complex exactly when the parity of the cardinality of the set of paths W in G verifying 
the following conditions is odd. Let W be a sequence of edges (ei, ...,efc). 

• All the edges of W are between generators of the same two Maslov gradings. 

• The path length k is odd. 

• Each intermediate vertex on the path is adjacent to an edge of the path that disappears 
before it in the homology. 

Our program uses this combinatorial characterization of <9 s hort m a very straightforward way. 
It finds all possible paths in the long oval complex graph, by using a breadth first search, and it 
counts them. A simple caching of the generators explored is done to avoid exploring dead-ends 
many times. The point of this method is that the huge graph explored doesn't have to be stored 
at any time. It makes this method very efficient in terms of memory. Of course, the search is quite 
slow, but this method is not called often enough for speed to be a problem. 

5.4. Technical aspects of the programming. Our program was written in Pythor@. We used 
Idle and Eclipse @ with PyDev [j as editing environment. Our program would have been much 

^http:/ /www. python. org/ 
^http:/ /www. eclipse. org/ 
^http: / / pydev.sourceforge.net / 
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After the homotopy: 
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Figure 9. The effect of the disappearance of a bigon. The small circles represent 
generators (a part of the corresponding generator is drawn beside each circle). The 
bipartite graph represents (a small part) of a boundary map between two consecu- 
tive Maslov gradings. The generators containing a disappearing intersection point 
disappear. We see that what was a path of length three between two generators, 
before the homotopy, becomes an edge between those two generators afterwards. 

slower without the wonderful program called PsyccH of Armin Rigo. Psyco is a kind of Just-In- 
Time compiler for Python which can apparently greatly accelerate (from two to a hundred time 
faster) almost any Python program. It enables programs in Python, a high level language, to run, 
even for the most computation intensive tasks, with speed comparable to programs in a low level 
language like C. 
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